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We use considerations of energy balance and dissipation to derive a self-consistent 
version of the shear-transformation-zone (STZ) theory of plastic deformation in 
amorphous solids. The theory is generalized to include arbitrary spatial orienta- 
tions of STZs. Continuum equations for elasto-plastic material and their energy 
balance properties are discussed. 

I. INTRODUCTION 

Important progress has been made in a recent series of papers by Falk, Langer and myself 
on the shear-transformation-zone (STZ) theory of plastic deformation in amorphous solids 

HQ. 

In the first paper in this series we introduced and explored an energetic approach 
to the STZ theory at temperatures far below glass transition temperature, which helped 
us to define the limits of the theory's form. The finite-temperature version of the theory 
developed in the second paper |2j allowed us to make predictions that were comparable to 
to experimental observations of the behavior of bulk metallic glasses (Kato et al.jsj], Lu et 
al.Q). The success and the questions that these studies posed prompt us to look more 
carefully at the fundamentals of the theory and understand the extent to which the simple 
approximations that we used were correct, and how to construct the theory without them. 
This paper is focused on further generalizing and expanding the low-temperature STZ theory 
of plasticity. In particular, we reexamine the physical significance of two parameters that 
occurred in the energy balance equations introduced in Ref. [lj; and we show explicitly 
how to derive the tensorial version of the theory, already used in Ref. §|, that is needed in 
order to describe situations in which the orientation of the stress changes as a function of 
position and time. Finally, for completeness, we derive a full set of elasto-plastic continuum 
equations of motion for this class of models. 

The STZ theory of plasticity of amorphous materials at low temperatures was proposed 



by Falk and Langer in [6|. It is based on the previous works of Cohen, Turnbull, Spaepen, 
Argon , which argued that non- crystalline solids plasticity is due to atomic rearrange- 

ments at localized sites. This picture has also been confirmed by a number of computational 
studies 00. However, unlike the earlier theories, the STZ theory focuses in detail on 
how rearrangements at the localized sites (shear-transformation-zones) occur, and identifies 
as important dynamical variables not only the concentration of the STZs, but also their 
orientations. This new variable allowed immediately to obtain a description of elastic and 
plastic behavior as an exchange of stability between the two steady states. Such a sim- 
ple mathematical treatment appears to us to be much more natural than the approach of 
traditional plasticity theory with its yield criteria. 

Moreover, the original STZ theory offered an explanation of a wide class of plasticity 
phenomena such as work hardening, strain softening, the Bauschinger effect, and others. 
But, as pointed out in 6], it had an inconsistency which implied that the proposed form 

n 

was not completely correct. The energetic approach introduced in flj allowed to correct the 
inconsistency for a simple case of quasilinear approximation. As shown there, even in such 
a simple form the STZ theory captured the important features observed in both mechanical 
tests and calorimetric measurements of glassy polymers at temperatures far below glass 
tuition texture 0. A generation of sneh an approach (also for the ^linear 
approximation) to higher temperatures [2j has proven to be quantitatively successful in the 
description of the viscoelastic response of bulk metallic glasses under tensile loading 
However, gued in 

urn 

, the application of quasilinear theory is limited. Most 
notably, the quasilinear approximation exaggerates plastic flow at small stresses and low 
temperatures, and reduces memory effects. As only the non-linear STZ theory can be 
expected to adequately describe molecular rearrangements, it must be further developed in 
order to reach precise quantitative agreement with experiment. One of the purposes of this 
paper is to expand the energetic approach introduced in [1] to the non-linear STZ theory. 
A major challenge in developing the STZ theory was defining the form of the STZ creation 

] 



and annihilation rates. In the original paper 



these rates were proposed to be proportional 



to the rate of plastic work ae p i. Since this work can become negative, it was obvious that this 
form was not acceptable (the inconsistency noted above). An easy (but artificial) remedy 
was proposed in p - to make them proportional to the absolute value |cx£p;|. This was 
sufficient for handling complicated numerical simulations of necking where ae p i becomes 
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negative during unloading. These simulations also explicitly demonstrated that thinning of 
the neck could continue even after stretching of the sample had been stopped. This raised a 
question whether the proposed forms of the STZ theory agreed with fundamental physical 
principles - the first and second laws of thermodynamics. It appeared that, indeed, they 
did - the plastic deformation of the neck was driven by the energy stored in the bulk of 
material, and this process was dissipative. 

Beyond that, the involvement of energy concepts in the consideration of the theory opened 
a different perspective. In this paper we make a conjecture that will be the basis of all of 
the following discussion - that creation and annihilation rates are proportional to the rate 
of energy dissipation. This conjecture of proportionality allows us to self-consistently define 
all components of the theory (Section The formalism developed there is a useful tool in 
limiting the arbitrariness of possible forms of the dynamical equations, transcending the cur- 
rent framework of low-temperature STZ theory. In Section IIHI we demonstrate conclusions 
of Section [H] on two important examples. 

Another significant limitation of the original STZ theory was that it considered STZs 
oriented in a single direction only. In earlier studies that had to deal with stress changing 
its direction the form of the theory for amorphous material, isotropic in its 

nature, had to be guessed on a phenomenological basis. In Section llVl of this paper we return 
to the microscopic basics and construct a theory that includes STZs oriented in all possible 
directions. Thereafter, we introduce an approximation that allows us to rewrite the theory 
in a simpler tensorial form, with order parameters being the first and second moments of the 
orientational density of the STZs. This tensorial form is comparable to the above mentioned 
phenomenological theory. 

In Section we combine ideas of the previous sections, applying the energetic approach 
from Section |n] to the isotropic model of STZ theory from Section HVl 

An understanding of energetic processes in the plastic degrees of freedom allows us to 
deal more carefully with spatially distributed systems, which we discuss in Section ED Here 
we put all of the ingredients together and write dynamical equations for an elasto-plastic 
material in two dimensions, preserving a clear picture of energy balance. 

In Section I VI II we present some arguments in favor of our conjecture of proportionality 
between the rate of creation and annihilation of STZs and the rate of energy dissipation. 
We also discuss some details that have been left out so far, but still may be important to 
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obtain quantitative agreement with experiment. 

II. ENERGY CONCEPTS IN THE STZ THEORY OF PLASTICITY 

The basic premise of the STZ theory is that the process of plastic deformation in an 
amorphous material is due to non-affine rearrangements of its particles in certain regions, 
that are called shear transformation zones. The original STZ theory simplistically considered 
all STZs as oriented in a single preferred direction. A two-dimensional sample was subjected 
to pure shear loading with a principal axis of the deviatoric stress tensor coinciding with the 
preferred direction. Throughout this section we will adhere to the same propositions. 

To be specific, we will call the zones elongated along the y-axis as "+" zones and the 
zones elongated along the x-axis as "— " zones. We will denote the density of zones in the 
"+" state by n+, and in the "— " state by n_. For pure shear the deviatoric stress tensor 
has the form: s xx = —s, s yy = s, s xy = 0. 

Following we can think of the plastic strain rate as the result of transitions between 
the states of STZs: 

Epi = Xv (R_n_ — R+n + ) , (II. 1) 

where e p i is the ^-component of the plastic strain rate tensor, R + is the rate of transitions 
from "+" to "— " states, R- is the rate of transitions from "— " to "+" states, A is the 
elementary increment of the shear strain, and v is a volume of the order of the STZ volume. 
Generally transition rates are functions of stress s or, equivalently, of the dimensionless 
variable s/ft, where ft can be interpreted as a sensitivity modulus[f|. This modulus has 
dimension of stress or energy density. Equation (jll.lj) also implies that all STZs have the 
same size, and therefore the constants A and v are the same for all zones. 

We suppose that STZs can also be annihilated and created, with the annihilation rate 
R a and creation rate R c . The creation rate, unlike transition and annihilation rates, can be 
understood only as a quantity defined per unit volume. Thus, we have: 

n± = R T n T — R±n± — R a n± + R c . (H-2) 

We can rewrite Eqs. (jll.lj) . ()H.2j) in a more convenient form. If we introduce a parameter 
To that specifies some time scale for transitions, and define rate functions S = t (R_ — R + )/2, 
C = r (i?_ + R+)/2, T = S/C, T = r Q R a , densities = 2R C /R a , n tot = n + + n_, 
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n& = n + — n_ } and dimensionless quantity e = Xvu^, we get: 

2/100X0 . -p /TT Q \ 

rohtot = r(noo - n tot ) , (II.4) 

roipi = — (Tntot - n A ) ■ (H-5) 

This system of equations is completely determined if we define the functions R±, R a and R c , 
which was done in (|. In this paper we will postpone choosing specific forms of the transition 
rates and corresponding functions C, S, and first focus on the creation and annihilation rates. 

From ()H.2j) we see that an important feature of this theory is that creation and anni- 
hilation of STZs are independent of their orientations and occur with equal probability for 
both orientations. This is not a completely trivial assumption. We disregard the possibility 
that creation and particularly annihilation can happen in connection with transition pro- 
cesses, and thus be more intense for one orientation of STZs than the other. However, the 
assumption that the creation and annihilation rates are independent of orientation is simple 
and plausible. Another observation we can make is that the creation rate is very likely to 
depend on the structure of material, or in other words, on such characteristics as packing 
fraction, free volume or structural disorder, as this rate is not only a dynamical, but also a 
structural property. This is also expressed in the fact that we can define it per volume of 
material, but not per STZ. On the other hand, the annihilation rate, as well as the transition 
rates, is less likely to depend on the structure of material. This is expressed in the fact that 
they can be defined as rates per STZ, and can be thought of as properties of STZs, but 
not the surrounding material, the influence of which on individual STZs can be described 
by averaged quantities, such as average stress. In further discussion we will assume that 
changes in the structure of material can be described by changes in STZ degrees of freedom 
only. Thus, our only internal dynamical variables are n to t and tia, while is assumed to 
be a constant. 

It was proposed in [f| to make the rates of creation and annihilation proportional to the 
rate of plastic work 2se p i. A peculiarity of this expression mentioned earlier is that these 
rates, by definition always positive quantities, can become negative. This happens because 
plastic work does not entirely dissipate. 

In general, the rate of plastic work done on a system can be represented in the form 

2se pl = + Q , (II.6) 
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where ip is the energy that is stored in the plastic degrees of freedom and in principle can 
be recovered, and Q is the dissipation rate - a non-negative function of stresses and internal 
variables. 

As annihilation and creation rates themselves are non-negative, we propose to make them 
proportional to the rate of dissipation Q. We will give some reasons why this proportionality 
can be true in section IVII1 but at the moment this proposition should be viewed as a 
conjecture that provides a physically sensible model and adequately describes mechanical 
and thermodynamical phenomena in amorphous solids. 

Now we are in a position to derive formulas for Q, R a and R c . We write Q = AR a = 
AT/tq, where A is a coefficient determining the proportion in which dissipated energy drives 
creation and annihilation rates. Generally, this coefficient can be a function of total STZ 
density n tot , but not ua, meaning that dissipation produces creations and annihilations of 
STZs independently of their average orientation already present in the sample. Later we will 
refine our conjecture and postulate that the annihilation and creation rates are proportional 
to the rate of energy dissipation not simply per volume, but per STZ. Thus, the coefficient 
A will be proportional to n to t- As the energy ip depends only on the internal variables wa 
and n tot , we have: 

# dip . dip . 

-rr = jr-nA + w—ntot ■ (II- 7) 

at oua on tot 

Then, using (jll.3|) . (jll.4|) and (jll.5|) . we derive from (jll.6j) : 



A-nAH- + (noo-n tot )^ 



(II.8) 



'<9n A ' v 00 tot 'dntot 

In (|11.8|) we must choose ip in such a way that T is always non-negative. If we look at V 
as a function of s, we conclude that both the numerator and the denominator must always 
be positive independently. The numerator is guaranteed to be positive if its two factors 
always become zero simultaneously, that is at Sq = fiT~ 1 {nA/n tot ), where T is assumed to 
be monotonic, and T" 1 is the inverse function of T. This gives: 

dip eo/V-if / \ mn\ 
~ — = — T (nA/ntot) ■ (H.9) 

on a f^oo 

From (jll.9|) . it follows that ip as a function of tt-a is defined uniquely. If we suppose that the 
energy ip must be extensive in n to t, we get: 

^ = 600— , (11.10) 

noo V \ n totJ J 
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where P(£) = J T~ 1 (x) dx and k is a constant. The term proportional to k plays an 
interesting role here. It determines how much energy is stored in the material due to the 
presence of the STZs. This energy can be recovered if the sample is annealed and thus the 
number of STZs is reduced. However, in the low-temperature theory we do not have any 
way to reduce the density of STZs if it is less than (see Eq. (jll.4)) ). Therefore, if we are 
conducting mechanical tests only, this part of the energy appears to be dissipative, although 
in general it is not. 

Now we refine our conjecture and postulate that the annihilation and creation rates are 
proportional to the dissipation rate per STZ. We can rewrite our equations in a simpler form 
by defining A = rttat/n^ A = n A /n OQ , A = aeofmtot/n^, V = ip/e p,, § = s/ft. Equations 
fJlLSjl . dm, jm, (|TLSj) and (jTT"To|l then give: 

r A = 2C(s)(T(s)A — A) — TA , (11.11) 

r A = T(l-A) , (11.12) 

r e pl = e C(s) (AT(s) - A) , (11.13) 
r _ 2C(g)(AT(g)-A)(g-T- 1 (A/A)) 

M(A, A) ' 1 J 

* = A(P(A/A)) + k) , (11.15) 

where the denominator of V is 

M(A, A) = aA - jT-i + (1 - A) (p + k\ . (11.16) 

n q 

In earlier papers |1|,I2|, where we used the quasilinear approximation, we chose a = 1 and 
k = 1/2. But these parameters have a physical significance and we will later study how 
their choice influences the behavior of material. 

Let us now look at the locus of the equilibrium points A = in the s-A plane (see Fig. 
[Q. The importance of these points is due to the fact that they determine the two states of 
the system - jammed and flowing. The line s = T _1 (A/A) is the locus of jammed states, 
here e v i = 0. The other solution, s = T -1 (A/A) + M(A, A)/A, is the locus of flowing states, 
where e p i is non-zero. Note the role that A is playing here. Its equilibrium value is equal 
to one. Accordingly, the lines plotted for A ^ 1 are not true equilibrium branches. We will 
call them quasi-equilibrium, as they change when A relaxes to one. 

The jammed and flowing branches can intersect only at the point where M(A, A) = 0. 
The dissipation rate also diverges at this point. In general, the value of the variable m = A/A 
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FIG. 1: General s-A diagram. Here it is plotted for specific parameters of quasilinear model 
(Sec. IIII A|) . but the diagram's topology is the same for the general case. The thick solid lines 
show two steady states - jammed(l) and flowing(2). The arrows show regions where A increases 
or decreases, which is determined by the sign of A from Eq. (IITTlT) . The thick dashed line (3) 
shows the saturation value of A, when A = 1. The three thin dashed and dash-dotted lines(4-6) 
are the mA-hnes for fixed a and the three different values of k. The thin solid lines(7-10) show 
quasi-equilibrium branches for those three values of k at some initial value of A. 

at this point is a function of A; we will denote it as m^. Because of the divergence in T, 
the dynamics of Eq. is such that m is always less than m\. Thus, the value of 

determines the maximum number of STZs that may flip in one direction; we will call it the 
saturation point. 

The function depends on the parameters a and k. Let us look at how their choice 
influences function's behavior. From Eq. (|11.16|) we find that when A = 0, rri\ = m is 
the solution of the equation n = mQT~ 1 (m ) — P(m ), and when A = 1, = rri\ is the 
solution of the equation a = m\T~ x (m\). What happens if a = a K = moT _1 (mo), so that 
mi = m ? We can check that in this case M(A, A) vanishes for any A, if m = m , meaning 
that m-A = m . Thus, we can formulate an important property of Eq. (IH.16|) : for any k there 
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is an a = a K , such that rri\ is independent of A. The behavior of the function is also 
simple, if a differs from a K . We can prove that if a > a K , the function is monotonically 
increasing, and if a < a K , m\ is monotonically decreasing. 

To illustrate different choices of parameters a and k, in Fig. we show plots of A = Km\ 
as functions of s = T _1 (mA), obtained by varying A, for fixed a and three different values 
of k. We will call such curves mA-lines; each of them is the locus of intersection points of 
the quasi-equilibrium jammed and flowing branches, when A varies. As the value of a is 
fixed, steady state branches coincide for different k when A = 1. Of the three values of k, 
the intermediate value is such that a K is equal to the given value of a. 

In an elasto-plastic material the total strain rate is given by 



where \i is the shear modulus (see also Eq. (|VI.9J) and the discussion thereof). Let us consider 
solutions of the system 11111. 13JII.17J) at a constant strain rate i to t- If the strain rate is 
small, the mA-lines coincide with the dynamical trajectories in the regime when A is evolving 
from some initial value A towards unity. In other words, the dynamical trajectory in the 
s-A plane first moves along the quasi-equilibrium jammed branch calculated for A = A (line 
7) until the intersection with the mA-line and then moves along the mA-line (for example, 
along line 4 for the smallest k). For higher strain rates the dynamical trajectories tend to 
lie to the right of the quasi-equilibrium jammed branch and the corresponding mA-line and 
evolve from zero to some point on the flowing branch at A = 1, determined by the value of 
itot- The final value of s can be smaller than intermediate values, thus producing a stress 
overshoot. As we can see from Fig. ^ a stress overshoot is more likely to happen for large 
values of k. For small values of k the stress increase is usually monotonic. 

It is hard to find compelling reasons why in a glassy material the saturation value m\ 
should be dependent on A. Therefore, we will suppose that for glasses a = a K . Indeed, this 
assumption produces behavior typical for glasses, such as essential strain rate dependence 
and stress overshoot. Further, we will also study a case when rri\ is dependent on A. 
This case may be relevant for description of polycrystals, clays or soils, if deformation in 
such systems is due mostly to rearrangement of individual crystals or grains, rather than 
deformation of grains themselves. The difference from glasses, to which we particularly wish 
to refer here, is the presence of an additional means of energy dissipation due to friction 



Slot = s P i + s/2fi , 
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between constituent particles, which we will model with a larger dissipation coefficient, that 
is, with a > a K . 

Now, to make the discussion clear and to put our previous works into the current more 
general framework, we consider the simple case of what we call the quasilinear version of the 
STZ theory. Such an analysis was presented in much detail in albeit only with a = a K . 



III. EXAMPLES 

A. Quasilinear theory 

In the quasilinear theory the transition rate functions are supposed to be linear functions 
of the shear stress s. Namely, we assume that C(s) = 1, S(s) = s, so that T(s) = s, 
T-\0 = £, P(0 = £72- From jEQSJ we get 

^ = A(m 2 /2 + K), (III.l) 

where m = A/ A. The expression 16)1 for M becomes 

M(A, to) = A(a - k) + k - (A + l)m 2 /2. (III.2) 

We find that m\ = 2(a — k + (2k — a)/(A + 1)) and a K = 2k. 

In Q], we chose a = 2k = 1, so that |toa| = 1. Then Eq. (|II.14J) becomes: 



4A(AS - A) 



2 



(III.3) 



(1 + A)(A 2 - A 2 ) ' 

Using Eq. (jlll.3|) in the dynamic equations (jll.llllll.T3|) we find that non-flowing steady states 
occur at s = A/A < 1 and flowing steady states at s = (1 + A)/(2A) - (1 - A)A/(2A 2 ) > 1. 
The exchange of stability occurs at s — 1. This value can be naturally associated with the 
yield stress. 

We solve Eqs. (jII.llHH.T3l 111.17)) numerically at the constant strain rate and show the 
results in Fig. [2] (a, b). We plot s- A trajectories and stress-strain curves for three initial val- 
ues of A. When the initial number of STZs is small - the sample is annealed - a pronounced 
stress overshoot is observed. For quenched samples,that is, when the initial value of A is 
large, the stress overshoot disappears. As shown in |l(, the constant strain rate simulations 



of this model are qualitatively similar to the available experimental data 
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FIG. 2: Quasilinear model. s-A and stress-strain diagrams for two cases: m\ is constant (Fig. 
a, b), and m\ is monotonically increasing (Fig. c, d). Lines 1, 2 are the jammed and flowing 
steady-state branches at A = 1, line 3 is the rrtA-line. Figs, a, c show dynamical trajectories for 
the stress-strain curves from Figs, b, d respectively. 

Now we shall consider other choices of constants a and k. One artificial difficulty with 
the quasilinear approximation is that some choices of these constants lead to m\ larger than 
unity, thus allowing m to assume non-physical values. To satisfy the condition \m\ < 1, 
additional conditions must be imposed on the acceptable values of a and k. There are two 
regions for parameters a and k, where \m\ < 1 for all A : (1) when k < 1/2 and k < a < 2k, 
here a < a K , and (2) when a — k < 1/2 and 2k < a, here a > a K . 

In Fig. I21(c, d) we plot the results of simulation for a = 1/2, k = 1/20, that is, when they 
are in the second region. The strain rate is small, thus the steady state on the s-A diagram 
almost coincides with the intersection of two steady state branches, and, as discussed in 
Sec. |nj the dynamical trajectory follows along the quasi-equilibrium jammed branch and 
then along the mA-line. The stress-strain curve exhibits strain hardening as observed in 
polycrystals, soils or clays [17] . Because such a strain-rate curve exists in the limit of an 
infinitely small strain rate, it can be rate independent for many decades on the logarithmic 
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scale. During strain hardening almost all the energy goes to creating more STZs. As we 
already noted, we can not get this energy back in mechanical tests, so in this sense, such a 
regime can be considered to be dissipative. 

B. Non-linear STZ model 

The quasilinear model is very useful as a toy model because of its simplicity It allows us to 
proceed much further in analytical and, often, numerical calculations. But it is mainly useful 
only to gain qualitative insight into the underlying dynamics, not to look for quantitative 
predictions. The most important drawback of the quasilinear model is its exaggeration 
of plasticity at small stresses. This drawback can be traced to the form of function C(s) 
which is constant in the quasilinear approximation, but in reality is vanishingly small at 
small stresses. This property is also responsible for suppressing the dynamics of A at small 
stresses and, thus, for memory effects. 

The general derivation of section |H] suggests that the energetic approach to the fully 
non-linear model will give qualitatively the same results as those of the quasilinear model, 
while fixing inaccuracies of the latter. In this subsection we briefly illustrate this point. 

In a full STZ model C and S can be arbitrary non-linear functions of shear stress. The 
important fact to note is that, unlike what is assumed in the quasilinear approximation, the 
function |T(s)| = \(R- — R+)/(R+ + R-)\ is always less than unity and asymptotically ap- 
proaches it when s — > oo. This causes the function \T~ 1 (m)\ to diverge when |m| approaches 
unity. Thus, when the denominator (jII.16|) vanishes, the value of \m\ = |A/A| is always less 
than unity. It can not exceed unity at any values of parameters a and k, but it also cannot 
be equal to one. Value of m equal to one corresponds to the complete saturation - the case 
when all STZs are oriented in one direction. It is puzzling that the non-linear theory does 
not allow this. But we will see in section that this is what must happen, if we take into 
consideration that in amorphous materials the STZs are oriented arbitrarily. 

Next, to proceed with numerical calculations we make a particular choice of functions 
R + and We will assume that functions R± have the form offered in jfj], that is R± = 
exp{— /?exp (±s)} /tq, where (3 = V*/vf, Vf is the average free volume and V* is of order 
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FIG. 3: Nonlinear model. s-A and stress-strain diagrams for two cases: toa is constant (Fig. a, b), 
and mA is monotonically increasing (Fig. c, d). Lines 1, 2 are the jammed and flowing steady-state 
branches at A = 1, line 3 is the TOA-lfne. Figs, a, c show dynamical trajectories for the stress-strain 
curves from Figs, b, d respectively. 

of the average molecular volume. So we find that 

C{s) = exp (— (3 cosh s) cosh ((3 sinh s) , 
S(s) = exp (— /3 cosh s) sinh ((3 sinh s) , 

T(s) = tanh ((3 sinh s) . (III.4) 

Now we can find that the function T _1 (m) = arcsinh(arctanh(m)//3). It diverges loga- 
rithmically at m — ► ±1, but the function P(m) and consequently the plastic energy \l/ given 
by (jII.15|) are finite at |m| = 1. 

For the numerical simulations at a constant strain rate loading we used the parameter 
(3 = 6. In Fig. |21 (a, b) we have rri\ = 0.992 = const, which is obtained with a = 0.444, 
K = 1/3. The stress rate is large, so that the equilibrium point on the flowing branch (2) 
is far from the intersection of the jammed (1) and the flowing (2) lines. As in Fig. El ( a > 
c), we demonstrate the results for three different values of A . In Fig. El (c, d), we chose 
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a = 0.444, k = 1/15, so that the mA-line is monotonically increasing. The strain rate is 
chosen to be small so the steady state point almost coincides with the intersection of the 
jammed and flowing branches. Note that analytically we can calculate very little in the fully 
non-linear model, and even the numerical solution requires not simply solving the system 
of differential equations, but also numerically calculating the integral P(m) at every step. 
But the analysis from Sec. |H] predicts much of the solution's behavior just from knowing 
how the function behaves based on the numerical values of a and k. We see that the 
plots in Fig. |3] are very similar to the plots in Fig. El for the quasilinear model, since the 
s-A diagrams are topologically the same. 



IV. ISOTROPIC STZ MODEL OF PLASTICITY 



In this section we generalize the STZ model of plasticity to the case of arbitrary spatial 
orientations of STZs and arbitrary orientations of the stress. The first attempt to make such 
a generalization starting from microscopic basics was made by M. Falk (l4|, but was not 
quite complete. 

Here again we will consider a two-dimensional homogeneous sample under a pure shear. 
To be specific, we will classify STZs in relation to the direction of the x and y axes. We will 
specify that an STZ is in the "+" state if the angle between the direction of its elongation 
and the y-axis is smaller than it /4; when the same is true with respect to the x-axis, we 
will say that the STZ is in the "— " state (see Fig. HJ). Note, we suppose that the "+" and 
"— " orientations of the zone are perpendicular to each other. Deviations from a right angle 
should be described as fluctuations beyond the mean field theory, and therefore will not be 
considered here. We write the pure shear in the form s^- = sdfj, where is the direction of 
the principal axis of the stress tensor, s = a/ SijSij/2, and 

4 = 28*8}-^. (IV. 1) 

In the above equation ef is a unit vector in the direction 0, so that d% x = —d^ y = cos 20 and 
dty = dy X = sin 20. We measure the angle in the counterclockwise direction relative to 
the x-axis. For the purposes of this section we could have chosen the principal axes of the 
stress tensor to be oriented along the x and y axes, but as we will further want to generalize 
this discussion for the case of arbitrary temporal evolution of the stress, we suppose that 
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FIG. 4: Classification of STZs as being in "+" or "-" states. 

is arbitrary. 

Then we suppose that only the diagonal component of the shear stress tensor in the 
direction of the zone orientation (the projection of the shear stress tensor on that direction) 
is important for the dynamics of transitions between the states of this zone. Thus, for the 
dynamics of the STZ population we write: 

Xt = R-(Sa)Xa - R+(S«)X+ ~ RaXt + Rc , (IV.2) 
Xa = R +( S a)xt - R-(Sa)Xa ~ R aXa + R c , (IV.3) 

where xt * s the density of STZs in the "_)_"/"_" state oriented at an angle a relative to the 
x axis, and s a = sdf.~ a is the projection of the shear stress tensor on the direction a. At 
this moment the density xt * s defined for angles from — 7r/4 to 7r/4. Note that all STZs are 
included in this range due to the circular symmetry. 

Now we note that our classification of zones as "+" and "— " depends on the choice of the 
direction of x and y axes, which is arbitrary. If one zone is in the "+" state in relation to a 
particular direction, then it is in the "— " state in relation to the perpendicular direction, that 
is, Xa+x/2 = Xa- O ur dynamical equations should not depend on such arbitrariness; they 
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should give the same results independently of a reference direction. Therefore, Eq. (|IV.3|) 
for the angle a ± it /2 must be the same as Eq. (j!V.2|) for the angle a, and vice versa. Thus, 
we conclude that the following relation for transition rates must hold: 

R+(s a ) = R-(-s a ) . (IV.4) 

We suppose that transitions do not change the volume of material. Thus, we must 
describe the elementary change in strain by a traceless tensor, which in two dimensions is 
proportional to dfj. Again, we suppose that the magnitude of this elementary change is 
always the same, only its orientation can be different. In analogy with Section |H] we have: 

4 = \v r 4 d% {R + {s a ) X + a - R-(sa)Xa) da . (IV.5) 

J -IT /A 

The region of integration in (|IV.5j) is chosen to count every STZ only once. However, from 
the symmetry for the angles a and a + tt/2 we conclude that the same integral is correct 
with any limits of integration in the form [— tt/A + 7, 7r/4 + 7]. 

As in Section m we can introduce rate functions S = r (_R_ — R + )/2, C = r (R + + RJ)/2, 
T = S/C,T = r R a and also densities = xt + Xa» xt = xt ~ Xai Xoo = ™oo = 2i? c / R a , 
where, as earlier, tilde means stress rescaled by /2, that is s a = s a /fi, §ij = s^/p,, I = s/Jjl. 
These functions also have symmetry properties: S(s a ) = — 5(s Q+7r / 2 ), C{s a ) = C(s a +n/2), 
T(s a ) = -T(s a+7r/2 ), and x^ 1 = aftr/ 2 » xt = -Xa +7T / 2 - Using these variables we can 
rewrite Eqs. (JE2}, dES)) and (ITV3|l as: 

ro^g = — ^ 4 d%C{s a ) (r(5 a )xS?* - Xa) da , (IV.6) 

Xoo J—tt/4 

T 0X t = 2C{s a ) (T(s a )xf - Xa) - rx* , (iv.7) 

ToXa' = ^(Xoo ~ XT) ■ (IV.8) 

These equations are analogs of Eqs. (|II.3|) . (|II.4|) . ()II.5|) . but with arbitrary orientations of 
STZs. Instead of the number of STZs in two different states their variables are the densities 
of STZs with different orientations. As it is hard to deal with such equations, where Xa 
essentially plays the role of a distribution function, we further show that these equations 
can be simplified under sufficiently relaxed assumptions, and instead of the density Xa we 
can introduce its moments - scalar and tensor variables. 

Instead of the angular density we can introduce the total density of zones in a sample 
ritot = (2/7r) f^^x^da, the equation for which is easy to get by integrating Eq. (jIV.8|) . 
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Instead of xt we introduce the tensor = /^1 4 d^x^da. To get dynamical equations for 
we multiply Eq. (|IV.7|) by cfg and integrate it over a: 

ptt/4 



2 / d?.C(s a ) (T(s a ) X f - x £) da - . (IV.9) 

J -it/ 4 

An assumption we will make here is that initially Xq* does not depend on a. Then 
according to Eq. (jIV.8j) Xa* is independent of a at all later times. Next, in the integral 
(|IV.9|) we will approximate the function C(s a ) by a function C(s) that depends not on the 
projection of the shear stress tensor on a given direction, but on the principal value of the 
shear stress s. The only role that the function C played in the original paper [f| was to be 
responsible for memory effects. It was a vanishingly small function for small stresses, and 
thus effectively froze the internal variables in an unloaded sample, preserving information 
about the previous loading. Our approximation keeps such dynamics intact. Now the 
integral in Eq. (jIV.9|) can be calculated. Together with Eqs. (jIV.6|) . ()IV.8|1 our system 
becomes: 

r £fj = — C(t) (f(^n tot - nij ) , (IV.10) 
r n, = 2C(f) (m S fn tot - n i3 ) - Tn, , (IV.ll) 
Tohtot = r(noo - ntot) , (IV.12) 

where we denoted 

/it /A 
d6T{ls cos 2$) cos 28 . (IV. 13) 

-7r/4 

In the derivation of this system we used the previously discussed property that we can 
change the region of integration to any quadrant. In more detail the integrals in (jiV6j) and 
(|IV.9|) had been calculated as follows: 

»7r/4 ptt/A 



/7T/4 /'7r/4 
d"jT(s a )da = / dfjT(s cos 2(0 — a))da = 

/7r/4 /*7r/4 
df- e T^ cos 26)d6 = d% \ cos 26 T(f cos 20) dO , 
■7r/4 J-tt/4 



-7r/4 J-tt/4 

where dfj is equal to s^/s. 

Equations fjlV.lOjl . pV.llj) and (|IV.12jl give the description of plasticity in the isotropic 
generalization of the STZ theory. 
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V. THE PROPORTIONALITY HYPOTHESIS FOR THE ISOTROPIC STZ 

MODEL 

We now show how to expand the results of Section |H] for the isotropic case. Again, as we 
will want to generalize results of this section for the case of arbitrary temporal evolution of 
the stress, we suppose that the principal axes of the tensor Sjj do not necessarily coincide 
with the principal axes of the tensor n^. We write the plastic work done on a system as: 

■pi _ e o (n-z^ij \ _ dip(njj, n tot ) 



We will denote s xx = -s yy = s, s xy = s yx = r, n xx = -n yy = n A , n xy = n yx = ns, and 
the invariant of the tensor as n = (n\ + n 2 & ) 1 / 2 . The energy ip is now a function of three 
variables, so: 

dip _ dip dn A | di/j dn s | dip dn tot ^ ^ 

dt dn& dt dns dt dn to t dt 
As in Section HU we suppose that Q = aeofm to tT /Ton^. Writing (jV.ljl in components and 
then assuming that the energy ip can depend on n& and ris only through n, we find: 

r = (n tot Tm - n A ) (g ~ |^) + (n^f - n s ) (gr - gg) 

The rate function T = r -R a must always be positive. In analogy with Section m considering 
this expression as a function of stresses allows us to conclude that the numerator and the 
denominator of T must always be positive separately. For fixed tia, ns and n to t and varying 
s, t we want the numerator to pass through zero at a single point (sq, t ) and be positive 
elsewhere. The numerator becomes equal to zero when its first and third brackets are equal 
to zero. This happens for stresses s = n A s /(n to tT(s )) and r = riss / (n tot T(s )). Now 
we can express So and r as functions of the variables n A , ns and n to t only. Noting that 
s = (sq + To) 1//2 ' we fi 11 ^ that $o = A*T -1 (n/n to t). Substituting s in the expressions for s 
and r , we get s = Ji'T~ l {n/n tot )n A /n, r = JiT~ l {n/n tot )ns/n. 

If the second and the fourth brackets also pass through zero at this point, they will always 
have the same sign as the first and the third brackets correspondingly, ensuring positiveness 
of the numerator. Thus, from either the second or the fourth bracket we find: 
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Therefore, we find that the energy if) has the same form as in Section |H] Eq. (|H.10|) : 

if) = e -^(p(^- )+k) . (V.5) 



Uoo V \ n tot, 

Now we can write out our final result for the tensorial generalization of the low tempera- 
ture STZ theory of plasticity in the form analogous to Eqs. (jH.llHH.T6]) . If again we denote 
A = ritot/rioo, Aij = n^jn^, A = n/n^, we get: 

Tbeg = e C(l) (AT (iff - A,,) , (V.6) 
7-oAy = 2C(I) (aT(I)^ - A,,) - rAjj , (V.7) 

r A = T(l - A) , (V.8) 
C(l)(Af(l)^-/ a - - Ajj)(gjj - T'\A/A)A lJ /A) 

r " m(a7A) • (v ' 9) 

Expressions M(A, A) and \I/(A, A) = ip/e^Jx are the same as (jII.16|) . (jII.15|) with A replaced 
everywhere by A. 

Finally, we can compare our tensorial theory with arbitrary spatial orientations of STZs 
and arbitrary loading, derived in this section, with the limited STZ theory of Section ITD for 
STZs oriented only along two preferred axes and pure shear loading. If we consider pure 
shear in the generalized STZ theory of this section, we must assume that the principal axes 
of tensors s^- and Ay are the same. Thus, for pure shear Eqs. (|V.6HV.8|) become the same as 
Eqs. fill. 11IIII. T3j) . Therefore, the results of Section |H] hold for the STZ theory generalized 
here. We also note that in the discussion of Section |H] an important role was played by 
the saturation point - the value of m for which no further transitions were possible. In 
the isotropic case, when all STZs are switched in one direction, that is, when \x&\ = X toi ■• 
we can find from the expressions for n to t and of Sec. IIVI that n = n tot . For other 
orientational distributions, when < x tot a t least in some interval of angles, we have 
m = n/ritot < 1- However, the projection of the stress tensor on the directions in the narrow 
strips under angles ±7r/4 to the principal axes of the stress tensor is small, for any finite 
value of s, leading to \Xa \ < X tot & t least for those angles. Therefore, we must expect that 
the saturation point will be reached at m = n/n tot < 1, as has been assumed in Sec. IIII Bl 
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VI. CONTINUUM EQUATIONS AND ENERGY BALANCE 

The plasticity described by the STZ theory can be incorporated into a continuum theory 
that describes elastic and plastic behavior of viscoelastic solids using a general framework, 



discussed, for example, in 



plastic l 

HQ. 



A. The STZ theory of plasticity in a spatially inhomogeneous situation 

We start with the generalization of the isotropic STZ model of plasticity for a spatially 
inhomogeneous situation. 

To make the physical picture clear, we will now discuss details omitted for simplicity in 
the previous sections. Let us consider a small region of material, much smaller than the 
size of the sample, but much larger than individual atoms and inter-atomic distances. This 
region contains many STZs of all possible orientations, but from a macroscopic point of 
view it is infinitesimally small and is identified by its coordinates only. Thus, we are on a 
mesoscopic scale. 

As this region contains many STZs, we consider the average effect of transitions between 
their states (which are changes in the positions of atoms on the microscopic level) on this 
region as a small part of the sample. From this point of view the transition between the 
states of an STZ gives rise to a change of strain at the point where this region is. 

Further we will describe the material by what is called the referential description jl^ . 
Namely, suppose that we are sitting in the material coordinate system and then at some 
time t we freeze our frame of reference and describe the evolution of the material during 
an infinitesimally small time interval in this frozen frame of reference. We can see that the 
discussion of Section HVl is correct even for an inhomogeneous situation in the material frame 
of reference, when the coordinate system not only moves with the particular small region of 
material, but also rotates with it. In the referential frame of reference we must exclude the 
effect of translational and rotational motion to make sure that we consider the same region 
of material under the same angle. 

Thus, instead of the time derivative of angle dependent quantities Xa 5 the dot 
in the expressions (|IV.2|) . (jIV.3|) and later must denote a complete co-rotational derivative: 
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where t>j and u are the translational velocity and the angular speed of our region. When 
deriving Eq. (jIV.9j) the integral df^x^da gives the tensorial co-rotational derivative 

Vriij driij driij . TT . 

-jjf = -Qj~ + Vk ~g^ + n ik w kj ~ w ik n kj , (VI.2) 

where = \/2{dvi/ dxj — dvj/dxi) denotes the spin tensor. This co-rotational derivative 
must be used in place of riy in (jIV.9|) and further. Correspondingly, instead of the time 
derivative h to t in the expression (jIV.12|) we get the total derivative 

dn tot dn tot dn tot ,,, TO , 

^r = ^r + ^' (VL3) 

as the rotational part integrates out. Finally, in the referential frame of reference the time 
derivative of the small strain tensor e?. is equal to the rate of deformation tensor Df-. 
Considering the above, the system f)IV. 10HIV. 12|) becomes: 

roDfj = e / e (po/p)C(I) (f(%) S -fn tot -n^ , (VI.4) 

r ^ = 2C(f) (f{%) S -fn tot - m>j - Truj , (VI.5) 

T °~lf~ = r ^°° _ ntot ^ ' (VI.6) 
In (jVI.4|) we took into account that the elementary strain, which is due to a transition 
between STZ states, can depend on the local density of material p (po denotes some reference 
density). We will discuss this point later. 



B. Continuum theory of elasto-plastic deformation 

Here we write out a complete set of equations needed to describe arbitrary elasto-plastic 
deformation of material. We also make an effort to demonstrate the energy balance prop- 
erties of our system of equations. This question is certainly not new for a system with 
constitutive relations in the rate form. However, we consider it important to show how 
plasticity described by the STZ theory can be incorporated into such a framework. 

First, we need to assert that our set of equations contains general equations which are 
true for any material - the conservation of mass and momentum balance equations: 

d P , dvi n t\n i\ 

Tt +p dx- = °> (VL7) 

pa * = ' (VL8) 
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where dj is the acceleration of material points, which in an inertial coordinate system is 
equal to dvi/dt, and is the true stress. 

We describe the material properties by a set of constitutive equations, which also includes 
equations for internal variables. To describe a viscoelastic solid, we additively decompose 
the total strain rate tensor D\f = l/2(dvi/dxj + dvj/dxi) as the sum of elastic and plastic 
parts, which is true under the assumption that elastic strain is small: 

Df n;] ■ Dfj . (VI.9) 

We would like to describe elastic behavior of the material simply by Hooke's law, but 
since here we are dealing with large deformations of solids and our equations are in the rate 
form, we need to take into account at least to some extent the dependence of the elastic 
properties of the material on its density. As we will see, this is dictated by the conservation 
of elastic energy. It is convenient to postulate that the equation of state of the material is 
defined by a function F K : 

p=-Kf K (p /p)F K (p /p) , (VI.10) 

such that ficix) = F' K (x), F K (1) = 0, fxiX) = 1- In the above equation p is the true pressure 
and po is the reference density of the material, which is convenient (but not necessary for 
further discussion) to assume to be the density of the material at zero pressure. The spherical 
part of the elastic response is fully described by this equation and, in fact, K here is the 
bulk modulus. Now we can introduce the conjugate stress and strain measures (the strain 
measure is given implicitly, by defining only the rate of deformation): 

V = vl !k{pqI p) ; Du = D ii f K (p /p)p /p . (VI. 11) 

Then, according to ()VI.9j) and (|VI.7J) . we can write the rate form of Hooke's law as 

OS - 4f , (V.,2) 

which coincides with the usual form in the case of small deformations. Similarly, for the 
deviatoric part of elastic response we have: 

1 Va dev 

mfyev = ± 1±_ Ni.YS) 

v l3> 2/i Vt ' V ; 

where the conjugate stress and strain measures are: 

st = vfr/fMp) ; m) dev = (Df J ) d ™fMp) Po /p . (vi.i4) 
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The conservation of mass equation (| VI. 7)1 . the momentum balance equation (|VI.8|) . the 
constitutive equations ()VI.4|) . ()VI.9jl . ()VL13)1 . the equation of state (jVI.10)) . and the equa- 
tions for dynamics of internal variables (jVI5)) . ()VI.6)1 constitute a full system of equations, 
which describe elasto-plastic behavior of a material. Those equations possess the property 

nnn 

of frame indifference |1 8L I19L 1201] . In particular, we used this system in a simplified form for 
simulations of necking 

The energy balance equation can be derived from the momentum balance equation. This 
derivation is very well known for the balance of energy in a volume fixed in space, but is less 
known for the case we are interested in here, when the balance of energy is considered in the 
volume of material. By multiplying Eq. (JVI.8)) by V{ and integrating over some arbitrary 
material volume V we get: 

f ( dv?pd? p d (pfj v f , \ f 

l IV) yd^ + j*Tt2K + j a Tt— + ) dV = J m v -"" dS ' • (VU5) 

The factor p/po in front of the total derivative plays an important role. Without it we 
would not be able to move differentiation over time in front of the integral. But as p/po is 
the Jacobian of the transformation from the coordinate system Xi(t) to the reference state 
Xj(0), we can first change the variable of integration to Xj(0), then put the time derivative 
in front of the integral (instead of a total derivative we will only be left with a derivative 
over time), and finally we can change variables of integration back to Xi(t). This is a purely 
mathematical procedure. It can be physically interpreted in the following way: instead of 
integrating over the time varying volume dV, we integrate over the conserved mass pdV . 
We get: 

Z-dev\2 

+ 

Po 




+ —^\dV= (VI. 16) 



P QdV . 



'(V) Po 

Above we also supposed that the plastic work can be expressed as 

P i P ( dip{n i:j ,n tot ) 



D ^^ [ ^_^ + a y (VU7) 

Equation () VI. 16)1 shows that energy in a particular volume of material consists of kinetic, 
elastic and plastic parts. It is changed by the work of external forces and it also dissipates 
due to plastic processes. 
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An important example relevant to above discussion is the Kirchoff stress tensor a. 



OjjPq/p, which is often used in engineering applications |2l| and standard en gine ering software 
22j. This stress tensor is conjugate to the rate of deformation tensor |23|. We get such 
a formulation, if we set Fk(x) = lnx, fx = l/ x > ft* = V x - This formulation assumes the 
following equation of state: p = K(p/p ) In p/po- 

Now we return to the assumption ()VL17[) . For this equation to be valid, the plastic rate 
of deformation tensor must be dependent on the density of material. This dependency has 
already been explicitly introduced in (|VI.4|) . At this point it is convenient to generalize our 
description of plasticity and also take into account the possible dependence of transition 
rates on the local density of material, which we include in the definition of the stress tensor 



Sij: 



% = *trfe(P0/p)P0/p ■ (VL18) 



Then the density of the rate of plastic work D^cr^ can be expressed as a product of p/p 
and a function of s^, ny, n tat , but not density. The equation (jVI.17|) then follows; we used 
it in the form (|V.l)) in connection with our hypothesis of proportionality of the annihilation 
and creation rates to the dissipation rate. 



VII. DISCUSSION 



Now that we have postulated that the rates of STZ creations and annihilations are pro- 
portional to the rate of energy dissipation and shown how to derive dynamical equations, 
we will proceed with discussing physical mechanisms that can underlie this hypothesis, and 
possible directions to further develop the STZ theory. 

The real microscopic picture of plastic deformation is far more complicated than what 
we describe in our model, where the properties of material are determined by the behavior 
of STZs only. At present, we can only tell if an STZ exists by observing localized atomic 
rearrangements - transitions from one STZ state to the other. But in principle an STZ is a 
spot where transition is potentially possible. However, since we judge about presence of an 
STZ only after the fact of transition, it is impossible to say whether an STZ was annihilated 
or created if we did not see its transition; or, even if we saw its transition, it is impossible 
to say at what point in time the STZ was created or annihilated. 

Thus, for example, it is much easier to understand the energetic properties related to a 
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transition that already happened. When atoms in an STZ rearrange, an additional stress 
field is created around the place of rearrangement. It is in this field that the plastic energy 
ip is stored, and this energy is in principle recoverable during a reverse transition. 

However, how do we understand the energetic processes related to the elusive events 
of STZ creations and annihilations? Annihilations are easy to imagine as impossibility of 
reverse transition after the initial transition or a series of transitions. In this case we can 
say that the STZ has annihilated and the energy stored in it has dissipated. Hence we can 
see a direct connection between the dissipated energy and annihilation. 

Let us look further. Any transition at low temperatures is a transition from a higher 
energy state to a lower energy state. This transition and creation of the stress field around the 
STZ is accompanied by dissipation of energy equal to the difference between the energy levels. 
This difference, before being absorbed by thermostat, can cause significant local increase of 
kinetic energy and additional atomic rearrangements which, along with transitions of other 
STZs, can lead to creations of new STZs and annihilations of existing ones. Thus, the energy 
dissipation is again related to creations and annihilations. 

Another important problem is to consider other essential degrees of freedom describing 
the structure of material. As we mentioned earlier, can be especially sensitive to them. 
In a theory for elevated temperatures, it is the increasing temperature dependence of 
that gives calorimetric characteristics of glass transition. An interesting way to introduce a 
variable describing disorder in the structure of material was offered in where 12^ was 
assumed to depend on that new variable instead of directly on the temperature. 

In the complex and not yet fully understood picture of the microscopic mechanisms 
underlying plastic deformation in amorphous solids, the conjecture of proportionality offered 
in this paper is the simplest of what can be suggested for STZ creation and annihilation 
rates, and it can be useful beyond the current framework of the low temperature theory. 

Acknowledgments 

This research was primarily supported by U.S. Department of Energy Grant No. DE- 
FG03-99ER45762. I particularly wish to thank Jim Langer for constant attention to this 
work and useful comments during preparation of this manuscript, and Lance East gate for 
many useful suggestions. I would also like to thank Craig Maloney, Anthony Foglia and 



26 



Anael Lemaitre for helpful discussions. 



[1] J. S. Langer and L. Pechenik, Phys. Rev. E 68, 061507 (2003). 

[2] M. L. Falk, J. S. Langer, and L. Pechenik (2003), accepted for publication in Phys. Rev. E. 

[3] H. Kato, Y. Kawamura, A. Inoue, and H. S. Chen, Applied Phys. Lett. 73, 3665 (1998). 

[4] J. Lu, G. Ravichandran, and W. L. Johnson, Acta Mater. 51, 3429 (2003). 

[5] L. O. Eastgate, J. S. Langer, and L. Pechenik, Phys. Rev. Lett. 90, 045506 (2003). 

[6] M. L. Falk and J. S. Langer, Phys. Rev. E 57, 7192 (1998). 

[7] D. Turnbull and M. Cohen, J. Chem. Phys. 52, 3038 (1970). 

[8] F. Spaepen, Acta Metall. 25, 407 (1977). 

[9] A. S. Argon, Acta Metall. Mater. 27, 47 (1979). 
[10] D. Srolovitz, K. Maeda, V. Vitek, and T. Egami, Phil. Mag. A 44, 847 (1981). 
[11] D. Deng, A. S. Argon, and S. Yip, Phil. Trans. Roy. Soc. Lond. A 329, 549 (1989). 
[12] O. A. Hasan and M. C. Boyce, Polymer 34, 5085 (1993). 
[13] M. L. Falk and J. S. Langer, MRS Bull. 25, 40 (2000). 
[14] M. L. Falk, Ph.D. thesis, UCSB (1998). 
[15] J. S. Langer, Phys. Rev. E 62, 1351 (2000). 
[16] J. S. Langer, Phys. Rev. E 64, 011504 (2001). 
[17] J.-G. Zhu and J.-H. Yin, Can. Geotech. J. 37, 1272 (2000). 

[18] L. E. Malvern, Introduction to the Mechanics of a Continuous Medium (Prentice-Hall, Inc., 
New Jersey, 1969). 

[19] J. Lubliner, Plasticity Theory (Macmillan Publishing Company, New York, 1990). 

[20] J. G. Oldbroyd, Proc. Roy. Soc. (Lond.) A200, 523 (1950). 

[21] R. M. McMeeking and J. R. Rice, J. Solids Structures 11, 601 (1975). 

[22] ABAQUS Theory Manual (Hibbit, Karlsson & Sorensen, Inc., 1998). 

[23] R. Hill, J. Mech. Phys. Sol. 7, 209 (1959). 

[24] J. S. Langer (2004), cond-mat/0405330. 



